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Abstract 



A comprehensive analysis of tadpole-improved SU(2) lattice gauge theory is 
made. Simulations are done on isotropic and anisotropic lattices, with and 
without improvement. Two tadpole renormalization schemes are employed, 
one using average plaquettes, the other using mean links in Landau gauge. 
Simulations are done with spatial lattice spacings a s in the range of about 0.1- 
0.4 fm. Results are presented for the static quark potential, the renormalized 
lattice anisotropy at/a s (where at is the "temporal" lattice spacing), and for 
the scalar and tensor glueball masses. Tadpole improvement significantly 
reduces discretization errors in the static quark potential and in the scalar 
glueball mass, and results in very little renormalization of the bare anisotropy 
that is input to the action. We also find that tadpole improvement using 
mean links in Landau gauge results in smaller discretization errors in the 
scalar glueball mass (as well as in the static quark potential), compared to 
when average plaquettes are used. The possibility is also raised that further 
improvement in the scalar glueball mass may result when the coefficients of the 
operators which correct for discretization errors in the action are computed 
beyond tree level. 
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I. INTRODUCTION 



Simulations of lattice Quantum Chromodynamics (QCD) have undergone significant 
changes in the past few years, with a host of new actions under investigation, and with 
the viewpoint emerging that simulations on coarse lattice can yield reliable results. A major 
impetus for these changes was the realization that the large radiative corrections which oc- 
cur in many quantities in lattice theories have a common origin, coming from cutoff effects 
due to tadpole diagrams specific to lattice actions 0. 

Tadpole diagrams in lattice theories are induced by the nonlinear connection between 
the lattice link variables U^{x) and the continuum gauge fields A^x), 

U^{x) = e iaaA ^ x) . (1) 

The cutoff dependence of ultraviolet divergent tadpole diagrams spoils naive power counting 
in the lattice spacing a. Higher dimension operators that are generated by lattice actions 
generally induce large radiative corrections, their contributions being suppressed by powers 
of a s (a), rather than of a. Fortunately, there is now considerable evidence that the effects 
of tadpole diagrams can largely be removed with a simple mean field renormalization of the 
links 

U,(x) - (2) 

Uq 

where an operator dominated by short-distance fluctuations is used to determine uq. (Al- 
ternative approaches to the design of improved lattice actions include the construction of 
fixed point actions @, and the use of nonperturbative renormalization conditions ||.) 

Tadpole improvement has helped to revitalize interest in the Symanzik improvement 
program |4[], with a number of new more complex actions currently under investigation, 
partly with the goal of doing precision simulations on coarse lattices (for a review see Ref. 
0). One of the earliest applications of tadpole improvement was in the development of the 
nonrelativistic QCD (NRQCD) action for heavy quarks [fj. Tadpole improvement is now 
widely used in large scale simulations of many actions |7j], and plays an important role in 
current efforts to extract continuum results from simulations on fine lattices. 

The current interest in simulations on very coarse lattices was stimulated by the first 
study of tadpole-improved gluonic actions ||. Simulations of the static quark potential, and 
of the spin-average charmonium spectrum in NRQCD, using an improved gluonic action on 
lattices with spacings as large as 0.4 fm, showed discretization errors of only a few percent. 

This has led to more recent efforts to further optimize tadpole-improved actions. In 
particular, there has been considerable work to determine the optimal choice of operator to 
use in defining the mean field renormalization factor uq P|-|I"3|]. Most previous simulations 
have used the fourth root of the average plaquette wo,p for tadpole improvement where, in 
SU(N) gauge theory (on isotropic lattices) 

/ 1 \ 1/4 

u °' p = \N ReTlU py ■ ( 3 ) 

However simulations of the static quark potential M, of the quarkonium spectrum in NRQCD 
m 12 , 13] , and of a relativistic fermion action [Tl[] have demonstrated that discretization 



errors are further reduced when the mean link in Landau gauge Uo,l is used, where 



2 



U ,L 



0. 



(4) 



There has also been a rapid evolution of more complex, highly improved actions (see e.g. 
Refs. |f), 11,11, IB]). In particular, there has been a revival of interest in actions defined on 
anisotropic lattices, where the lattice spacing in the "temporal" direction a t is kept smaller 
than the spatial spacing a s [0,0 . This enables much more efficient simulations of hadronic 
systems with large masses, for example, where the exponential suppression of the correlation 
function becomes prohibitive on lattices with large a t . This is especially relevant for glueball 
simulations, where discretization errors may be acceptably small even with a s as large as 
0.4 fm |fl6|-i~8[j, but where the correlation function becomes extremely noisy after only a few 
time steps if a comparable a t is used. 

An impediment to the use of anisotropic lattices has been the need to measure the 
renormalized anisotropy (at/a s ) p h y s in the simulation which, in the case of the Wilson gluon 
action , can differ appreciably from the bare anisotropy £ that is input to the action 



v u> s / bare 



(5) 



However it has recently been shown that tadpole improvement reduces the renormalization of 
£ to a few percent [nj,|T5[| , which is small enough to be neglected in many applications. Recent 
simulations of a variety of glueball states with a tadpole-improved action, on anisotropic 
lattices with coarse a s , have provided results that compete with much larger scale simulations 
on fine lattices 



17,18 



In this paper we present a comprehensive analysis of the tadpole-improved SU(2) lattice 
gauge theory. We do simulations on both isotropic and anisotropic actions, with and without 
improvement. Results are obtained with two tadpole renormalization schemes, one using 
average plaquettes, and the other using mean links in Landau gauge. We also compare with 
simulations of the Wilson action. Simulations are done with spatial spacings a s in the range 
of about 0.1-0.4 fm. Results are presented for the static quark potential, the renormalized 
anisotropy, and for the scalar and tensor glueball masses. Some of our results have been 
reviewed in preliminary form elsewhere p,0,|l|,|l8|] . A full account of this work is presented 
here for the first time f2]| . (Simulations of the tadpole- improved SU(2) gauge theory on 
isotropic lattices have also recently been reported in Refs. |2"2||2"3| .) 

We note that there is a long history of simulations of SU(2) gauge theory, motivated 
by the fact that this theory exhibits much of the physics found in SU(3) color, including 
linear confinement of static quarks and a rich glueball spectrum. Simulations of SU(2) color 
have been used to shed light on the physics of confinement, and to test new algorithms for 
use in more realistic simulations. The reduced computational cost of simulations in SU(2) 
continues to be exploited in, for example, recent large scale simulations of the static quark 
potential pif 



The work we present here on tadpole- improved SU(2) gluonic actions is not only of 
interest in reproducing the features of the SU(3) theory, but also serves to suggest new 
avenues for further development. Of particular interest are new results presented here for 
the SU(2) scalar glueball mass, which relate to a peculiar feature that has previously been 



observed in SU(3) coarse lattice simulations [|17]]. We find that the SU(2) scalar glueball mass 
first decreases as the lattice spacing a s is increased, reaching a minimum at a s « 0.3 fm; the 
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mass then gradually increases with a s . A similar "dip" in the SU(3) scalar glueball mass 
was reported in Ref. [|17]]. It has been conjectured [[L7[l8 that this behavior may be related 
to the presence of a critical endpoint in a line of phase transitions in gluonic actions that 
include an adjoint coupling |25| . 

On the other hand, we find that the depth of the dip in the SU(2) theory with plaquette 
tadpole improvement is about half of that found with the Wilson action. We also find that 
the depth of the dip is further reduced when mean link tadpole improvement is used (the 



SU(3) simulations in Ref. |L7j were done with average plaquette tadpoles). This suggests 



that the dip may be due, at least in part, to discretization errors that are more fully removed 
with mean link tadpoles. The tensor glueball mass exhibits relatively small scaling violations 
even with the Wilson action, which could indicate that the tensor has a larger size than the 
scalar. In this connection, we note that the scalar glueball has been found to be much less 



sensitive to finite volume effects than the tensor ||26|| . This supports the conjecture that the 
scalar glueball has a smaller size than the tensor, and hence is more sensitive to discretization 
corrections to the action. 

We also performed simulations after making small changes to the coefficients of the 
operators which correct for discretization errors in the action. These results raise the pos- 
sibility that the dip in the scalar glueball mass might be eliminated once the correct 0(a s ) 
renormalizations of the relevant operators in the action are included. 

The rest of this paper is organized as follows. In Sect. [Q] we present results from an 
improved SU(2) action on isotropic lattices. Two independent sets of simulations are done, 
using the two tadpole renormalization schemes discussed above (Eqs. (||) and ([|)). We 
simulate on lattices with spacings of about 0.25 fm and 0.40 fm. We compare the static 
quark potential with improvement to the results of a simulation using the Wilson action. 
We also present results for the scalar glueball mass. Details of our Wilson loop and glueball 
correlation functions are also discussed in Sect. ||. 

In Sect. [TTT] we present results from our simulations on anisotropic lattices. We run at 
several spatial lattice spacings a s in the range of about 0.1 fm to 0.4 fm, where the temporal 
spacing a t is kept near 0.1 fm. We run simulations for improved actions using the two tadpole 
renormalization schemes, and we also run on several lattices using the Wilson action. We 
compare bare and renormalized anisotropies, static quark potentials, and scalar and tensor 
glueball masses from these actions. We also illustrate the sensitivity of the scalar glueball 
mass to small changes in the coefficients of the operators which correct for discretization 
errors in the action. We briefly summarize our results and outline possible directions for 
further work in Sect. IT 



II. ISOTROPIC LATTICES 

The isotropic tadpole-improved action for SU(2) is identical in form to the SU(3) action 
], and at tree level is given by 



o _ _a \ "* j 5 Pjiv 1 Rjiv 1 Run [ /„x 

^,13«S 12 ug 12 ul /' [b) 



where P^ v is one half the trace of the lxl Wilson loop in the // x v plane, and is 
one half the trace of the 2x1 rectangle in the [i x v plane. The notation used in Eq. (^) 
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differs slightly from that used in Ref. ||, where a factor of 5/3mq was absorbed into the 



definition of (3. The notation used here follows that introduced in Refs. [p4|,|15| . The leading 
discretization errors in this action are of 0(a s a 2 ) and 0(a 4 ). 

We ran simulations at two lattice couplings using plaquette improvement, u = u P (Eq. 
(0)), and at one lattice coupling using mean link improvement, uq = u ,l (Eq. (f|)). For 
comparison, results were also obtained for a standard Wilson action on a coarse lattice. The 
parameters of these four lattices are given in Table |. In the following we will use f3p and 
@l t° denote the lattice couplings for the improved action with plaquette and mean link 
tadpoles respectively, and (3w for the Wilson action coupling. 

Results for the static quark potential V(R) on the lattices with comparable spacings 
are shown in Fig. [j]. The potentials were measured at integer separations R/a = 1-4. 
Measurements were also made using non-planar Wilson loops, corresponding to separations 
R/a = \/2, \/3, v^5, and v^8- Symmetric combinations of the shortest spatial paths con- 
necting two lattice points were used in the non-planar Wilson loop calculations. The data 
in Fig. m were obtained by looking for a plateau in the time-dependent effective potential 



aV{R,T) = - In 



W{R,T) 



W(R,T-1) 



(7) 



where W(R,T) denotes the Wilson loop. An estimate of the systematic error in the large 
T extrapolation is given by the statistical error in the second or third time step after the 
onset of the plateau. 

An iterative fuzzing procedure was used to construct Wilson loop (and glueball) operators 
with a very high degree of overlap with the lowest-lying state. Fuzzy link variables U- (x) 
at the nth step of the iteration were obtained from a linear combination of the link and 
surrounding staples from the previous step \27\ 



Ut\x) = Ut 1 \x) + ei: U^ n - 1 \x)ut 1 \x + j)U j {n - 1) \^ + ^ (8) 

where i and j are purely spatial indices, and where the links were normalized to WU = I 
after each iteration. Wilson loops (and glueball correlators, Eq. fll2| ) below) were constructed 
as usual by using the fuzzy link variables in place of the original links. The optimal number 
of iterations n and the parameter e typically varied from (n, e) = (10, 0.04) on lattices with 
a ~ 0.4 fm to (n, e) = (5, 0.4) on lattices with a ~ 0.1 fm. 

A defect of the action Eq. (^J) is that the gluon propagator acquires a high energy pole 
with a negative residue (at least in perturbation theory), due to terms in the action with 
two links in the time direction |14|,|15| . One consequence of these high energy "doublers" is 
that correlation functions exhibit exponential decay only at sufficiently large times. This is 
illustrated by an increase in the time-dependent potential V(R,T) at small times, as shown 
in Fig. |2|. This complicates the extraction of masses on coarse isotropic lattices, given the 
severe suppression of the correlation function at large times. This problem can be avoided by 
working on anisotropic lattices with sufficiently small a t , where elimination of 0{a\) errors 
in the action may be unnecessary, as described in the next section. 

The solid lines in Fig. [I] show the results of a fit of the on-axis data to a standard infrared 
parameterization of the continuum potential: 
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V &t (R) = *R-^ + c (9) 

Infrared fluctuations of a one-dimensional string give b = n/12 [PHI - However much better 
fits were obtained by leaving b as a free parameter. We typically found 6 ^ 0.1 on the 



coarsest lattices, with the string value emerging on the finest lattices analyzed in Sect. [Ill 



(see also Ref. fll?]). The lattice spacings a were extracted by matching the fit values for the 
string tension in lattice units to a physical value of y/a = 0.44 GeV. 

Discretization errors in the lattice action break Lorentz symmetry, which is clearly visible 
in the off-axis potentials || . These errors are dramatically reduced when the action is tadpole 
improved, as seen in Fig. [T]. As a quantitative measure of the improvement, we compare the 
potential measured in the simulation from non-planar Wilson loops with an interpolation to 
the on-axis data: 

AV(R) ee (10) 
uR 

Results for R = (a, a, a) are given in Table |I[ With tadpole improvement the error is only a 
few percent even on lattices with spacings as large as 0.4 fm, compared to an error of about 
30-40% for the Wilson action M. Discretization errors are further reduced when tadpole 
improvement is done with Uo,l, compared to when m 0j p is used. 
Glueball correlation functions G(T) were also calculated 

G(T) = <0(T)0(O)> - <0> 2 , (11) 

where 4>(T) denotes a linear combination of the trace of 1 x 1 (fuzzy) plaquettes Py(x, T), 
summed over spatial positions. We considered linear combinations which, in the continuum 
limit, excite scalar (J p = + ) and tensor (J p = 2 + ) glueballs: 

<Po+(T) = ]T [P 12 (f, T) + P 13 (x, T) + P 23 (x, T)] , 

X 

<MT) = £ [Pu(x,T) - P 13 (x,T)} , (12) 

X 

together with two other linearly independent combinations of plaquettes for the tensor chan- 
nel. Effective masses m(T) were computed in the usual way: 

An effective mass plot for the scalar glueball on the improved lattice with a ~ 0.26 fm 
is shown in Fig. |3|. The data are very noisy even at T of two or three lattice units, due to 
the large lattice spacing and glueball mass, despite the fact that an ensemble of over 30,000 
configurations was generated. The data for the tensor glueball were too noisy to be of use. 
As with the SU(3) simulations reported in Ref. |16[], the extraction of the glueball mass from 



a plateau in m(T) is computationally demanding on coarse isotropic lattices. As shown in 



the next section a far more efficient approach, proposed in Refs. is to simulate on 



anisotropic lattices with small spacings in the time direction. On the other hand, the results 



in Fig. |3| compare favorably with Wilson action calculations on fine lattices [29]. In Ref. 



| 30| a continuum extrapolation of Wilson action glueball data was made, with the result 



m Q +/y/a = 3.87 ±0.12. 
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III. ANISOTROPIC LATTICES 



The tadpole-improved SU(2) action on anisotropic lattices is again identical in form 
to the SU(3) theory and, following Refs. [I4|,15|, we omit corrections for 0(a 2 ) errors by 



working on lattices with small a t . This has the advantage of eliminating a negative residue 
high energy pole in the gluon propagator. The resulting action has rectangles R s t that 
extend only one lattice spacing in the time direction: 



S 



-P £ e{ 

x,s>s' \ 



5 P ss > 1 R ss > 1 R s t s 



3 u\ 12 u% 12 ul 
4 P st 1 Rst 



u 2 s u 2 t 12 u*v$ 



(14) 



where £ is the bare anisotropy, Eq. (^j). "Diagonal" correlation functions computed from 
this action decrease monotonically with time. 

On an anisotropic lattice one has two mean fields u t and u s . A natural way to determine 
them is to use the mean links in Landau gauge || 

u t = (iReTr U A ) , u s = (|ReTr U s ) , (15) 

where a lattice version of the continuum Landau gauge condition d^A^ = is obtained by 
maximizing the quantity 

£_L ReT r[/ M (x). (16) 

x,n u V a fj, 

Alternatively one can define the mean fields using the measured values of the average pla- 
quettes. One possibility is to first compute u s from spatial plaquettes, u A s = P ss i, and then 
to compute u t from temporal plaquettes, u 2 u 2 = P st . However, in the limit that a t /a s — ► 0, 
this procedure yields u t — > l/u s [^TJ; with mean link tadpoles Eq. ( |i5|) on the other hand 
one has the more physical limit u t = 1 — 0((a t /a s ) 2 ). Since the lattice spacing a t in our 
simulations is small, we adopt the following prescription p4]JT5|JT7[] for the mean fields in 
"plaquette improvement" 



U S = (Pss>) 1/4 . (17) 



In this connection, we note that the factor u t in the action Eq. (|14|) can in fact be absorbed 
into a redefinition of (3 and of the input anisotropy, according to (3/u t — > (3, and £,u t — > ^. 
As shown below one finds very little renormalization of the input anisotropy £ with either 
of the tadpole renormalization schemes Eqs. ( |15|) or (|i~7D. 

Simulations were performed on four lattices with plaquette improvement, and three lat- 
tices with mean link tadpoles. In addition, we ran simulations on six lattices with the Wilson 
action. The parameters of these thirteen lattices are given in Table [TT]. The spatial spacings 
a s lie in the range of about 0.1 fm to 0.4 fm, with temporal spacings a t kept near 0.1 fm. 

Configurations were generated using a heat bath algorithm. The number of updates 
between measurements varied from 10 on the coarsest lattices to 20 on the finest; integrated 
autocorrelation times were computed, and satisfied r int < 0.5 in all cases. Ensembles of 
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about 2,000 configurations were used to measure the static quark potential, while between 
80,000 and 160,000 configurations were generated at each f3 for glueball measurements. 

Tadpole improvement eliminates most of the renormalization of the input anisotropy £, 
as was first shown in Refs. [14],|l5||. This can be demonstrated by comparing the static quark 



potential computed from Wilson loops W x t, where the time axis is taken in the direction of 
small lattice spacings a t , 

W xt (R = ni a s ,T = n 2 a t ) nl^o Z xt e~ n ™ v{ - R \ (18) 

with the potential computed from Wilson loops W xy with both axes taken in the direction 
of large lattice spacings a s 

W xy (R = n 1 a s ,T = n 2 a s ) r^Zc Z xy e~ n2asV{R l (19) 

The physical anisotropy is determined after an unphysical constant is removed from the 
potentials, by subtraction of the simulation results at two different radii 

o*\ _ a t V xt (R 2 ) - a t V x t{Ri) ^ 

a sJ phys a s V x y(R2) — a s V xy (Ri) 

where a t V xt and a s V xy are the potentials in the lattice units relevant to the two sets of Wilson 
loops W xt and W xy , respectively (alternative methods for determining the renormalized 
anisotropy are analyzed in Ref. [pOfl). 

We computed the anisotropy twice, using two different radii R\ for the subtraction, with 
fixed i?2 — 2a s . We compare the anisotropies determined with Ri = a and \/2a in Table [II. 



Although it is advantageous to use the potential at smaller R, where the statistical errors are 
smaller, it is important to check that we are not sensitive to possible discretization errors of 
0(a^/R 4 ) [ftlj] . The two determinations of the anisotropy are in fact in excellent agreement. 
These results show that the input anisotropy £ is renormalized by less than a few percent 



when the action is tadpole improved |14| , |15|1 , at least over the wide range of lattices analyzed 
here. The renormalization is small enough that one need not measure the anisotropy in 
many applications. This is to be contrasted with the Wilson action, where the measured 
value of a t /a s is found to be about 20% lower than £ on the lattices analyzed here. 

Sample results for the potentials both with and without tadpole improvement are shown 
in Fig. The renormalization of the input anisotropy in the case of the Wilson action 
is plainly visible as a difference in slope of the potentials computed from W xt and W xy . 
A significant reduction in rotational symmetry breaking with the improved action is also 
apparent in the off-axis potentials, as illustrated in Fig. 

The scalar and tensor glueball correlators were calculated on all thirteen lattices. Repre- 
sentative effective mass plots are shown in Fig. |5|. The fuzzy correlators described in Sect. |Tl 



(Eqs. (H) and fll2l) ) were used. We note that the efficiency of these calculations could be 



improved by using a variational basis of several operators, rather than just the fuzzy lxl 
plaquette used here [fT7 |. Our final results for the masses were obtained from single expo- 
nential fits to the correlation functions. In most cases acceptable fit results for the scalar 
glueball were obtained from the fit interval T = (3-5)a t (we required that these fits be in 
good agreement with fits using larger intervals in T). The results are given in Table |TV| . 

The scalar and tensor glueball masses are shown as functions of lattice spacing squared 
in Figs. |6| and [5]. The tensor glueball data exhibit relatively small scaling violations, even 
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for the Wilson action on the coarsest lattices studied here. This could indicate that the 
tensor glueball has a very large size. It would be interesting to compare with the situation 
for the SU(3) tensor glueball, but Wilson action data are not available on sufficiently coarse 
lattices (for a compilation see Ref. |T7|). 

The SU(2) scalar glueball data exhibit a peculiar feature that has previously been ob- 
served in SU(3) coarse lattice simulations |17]]. The scalar glueball mass first decreases as 
a s increases, reaching a minimum at a s « 0.3 fm; the mass then gradually increases with a s . 
It has been conjectured [p/^ ,l8| that this "dip" may be related to the presence of a critical 



endpoint in a line of phase transitions in gluonic actions that include an adjoint coupling 



The data obtained here shed new light on this behavior. In particular, we find that 
the depth of the dip is reduced by about half when the plaquette-improved action is used, 
compared to results obtained with the Wilson action. We also find that the depth of the 
dip is further reduced when mean link improvement is used. This suggests that the dip may 
be due, at least in part, to discretization errors that are more fully removed with mean link 
tadpoles. 

We note that the SU(3) simulations in Ref. Jl7j were done with plaquette improvement. 
The SU(3) scalar glueball dip is more pronounced than in SU(2) with plaquette improvement, 
but Wilson action data are not available on sufficiently coarse lattices to allow a trend to 
be established in the SU(3) theory. Simulations with mean link improvement for SU(3) 
glueballs have not been done. One the other hand, some SU(3) simulations have been done 
with an adhoc change to the action, designed to move the theory away from the critical 
endpoint in the fundamental- adjoint plane |18[]; these results suggest that this mechanism 



does play a role in the SU(3) scalar glueball dip. 

In this connection we performed one last set of simulations, after making an adhoc change 
to the coefficient of the rectangle terms R ss > and R s t in Eq. (|T4]). We multiplied these terms 
by a factor of 1.2, which is consistent with an 0(a s ) renormalization of these operators on 
lattices with spacings in the range considered here. We ran at (3l = 0.85, with £ = 0.333 
(after retuning the Landau gauge mean link factors, where u t = 0.975 and u s = 0.758). The 
lattice spacing a s = 0.314(1) fm is in the middle of the dip found with the other actions. 
We find glueball masses m + / \fa = 3.82(7), and m 2 + / \/a = 5.5(4) (where a measured 
renormalization of £ by about 8% is included). These results are shown as an open triangle 
in Fig. H and [7|. These results raise the possibility that the dip in the scalar glueball mass 
might be eliminated once the correct 0(a s ) renormalizations of the operators which correct 
for discretization errors in the action are included. This conjecture is supported by finite 
volume studies of the glueball spectrum, which show that the scalar is much less sensitive 



to finite volume effects than the tensor p6| 



IV. SUMMARY AND OUTLOOK 



Tadpole improved SU(2) lattice gauge theory was applied to calculations of the heavy 
quark potential, the renormalized lattice anisotropy, and the scalar and tensor glueball 
masses. We analyzed improved actions on isotropic and anisotropic lattices, and we com- 
pared simulations with mean link and plaquette tadpole improvement. Comparisons were 
also made with simulations of the Wilson action. Tadpole improvement significantly reduces 
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discretization errors in the static quark potential, and results in very little renormalization of 
the input anisotropy. We also found a "dip" structure in the scalar glueball mass, analyzed 
as a function of lattice spacing, which has previously been observed in SU(3) simulations. 
We found evidence that this dip may be related to discretization errors in the action that 
are more fully removed when tadpole improvement is done using the mean link in Landau 
gauge. The possibility was also raised that further improvement in the scalar glueball mass 
may result when the coefficients of the operators which correct for discretization errors in 
the action are computed beyond tree level. 

Simulations with mean link improvement (as well as of the Wilson action) on coarse 
lattices in the SU(3) theory could help to clarify the origin of the dip in the scalar glueball 
mass. It might also be fruitful to analyze an improved SU(2) theory with an adjoint cou- 
pling included. A calculation of the scalar and tensor glueball sizes would also yield useful 
information, with the reduced computational cost of simulations in the SU(2) gauge theory 
providing further incentive for continued study of this system. 
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TABLES 



Action 


P 


(iReTrC/ pl )V4 


(iReTrig 1 /* 


a (fm) 


Volume 


AV{V3a) 


Improved 


0.730 


0.844 


0.779 


0.394(2) 


8 4 


0.04(1) 


(Plaquette) 


0.935 


0.872 


0.817 


0.264(1) 


8 4 


0.02(1) 


Improved 


0.550 


0.845 


0.779 


0.404(2) 


8 4 


0.01(1) 


(Landau) 














Wilson 


1.700 


0.802 


0.742 


0.402(2) 


8 4 


0.32(4) 



TABLE I. Simulation parameters for three isotropic tadpole-improved lattices, and one Wilson 
action lattice. The mean field for plaquette improvement (second column) and mean link improve- 
ment (third column) are given. The lattice spacings were determined from the string tension, with 
y/a = 0.44 GeV. The measured errors AV in the off-axis potential at R = (a, a, a) are also shown 
(see Eq. @). 



Action 


& 




u t 


u s 


a (fm) 


Volume 


Improved 


0.848 


0.276 


1. 


0.793 


0.366(1) 


8 3 x 32 


(Plaquette) 


0.981 


0.333 


1. 


0.819 


0.298(1) 


8 3 x 24 




1.114 


0.409 


1. 


0.842 


0.238(1) 


8 3 x 20 




1.214 


0.500 


1. 


0.860 


0.202(2) 


10 3 x 20 


Improved 


0.650 


0.276 


0.977 


0.718 


0.384(2) 


8 3 x 32 


(Landau) 


0.795 


0.333 


0.974 


0.757 


0.304(1) 


8 3 x 24 




0.905 


0.409 


0.966 


0.785 


0.244(1) 


8 3 x 20 


Wilson 


1.950 


0.250 






0.393(2) 


8 3 x 32 




2.000 


0.250 






0.355(2) 


8 3 x 32 




2.140 


0.333 






0.308(2) 


8 3 x 24 




2.243 


0.400 






0.232(2) 


8 3 x 20 




2.300 


0.500 






0.203(1) 


10 3 x 20 




2.400 


1. 






0.128(2) 


12 4 



TABLE II. Simulation parameters for the lattices analyzed in Sect. [II. The bare anisotropies 
£ and the mean fields ut and u s for tadpole improvement are shown, along with the spatial lattice 
spacings a s determined from the string tension. 
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{&tl Gs)phys 



Action 


P 




(R x = a s ) 


(R x = V2a s 


) AV(V3a) 


Improved 


0.848 


0.276 


0.273(4) 


0.272(8) 


0.08(1) 


(Plaquette) 


0.981 


0.333 


0.335(5) 


0.332(10) 


0.05(1) 




1.114 


0.409 


0.413(3) 


0.418(8) 


0.03(1) 




1.214 


0.500 


0.502(3) 


0.506(8) 


0.03(1) 


Improved 


0.650 


0.276 


0.277(11) 


0.269(22) 


0.06(1) 


(Landau) 


0.795 


0.333 


0.335(4) 


0.330(9) 


0.04(1) 




0.905 


0.409 


0.408(8) 


0.402(12) 


0.03(1) 


Wilson 


1.950 


0.250 


0.189(2) 


0.183(7) 


0.19(1) 




2.000 


0.250 


0.200(5) 


0.214(13) 


0.15(1) 




2.140 


0.333 


0.267(4) 


0.271(7) 


0.11(1) 




2.243 


0.400 


0.338(6) 


0.335(16) 


0.07(1) 




2.300 


0.500 


0.428(3) 


0.419(9) 


0.06(1) 




2.400 


1. 






0.04(1) 


TABLE III. Measured anisotropies (at/a s ) p ky S compared to the input anisotropies £ for the 


three actions. Each anisotropy was measured twice, using two different radii R\ 


for the subtraction, 


with fixed R s = 2a 


s (see Eq. © 


). The 


errors AV in the off-axis 


potential are also given. 


Action 







a s (fm) 




m 2 +l\[o 


Improved 


0.848 




0.366(1) 


3.49(7) 


4.4(7) 


(Plaquette) 


0.981 




0.298(1) 


3.41(6) 


5.6(3) 




1.114 




0.238(1) 


3.55(5) 


5.9(1) 




1.214 




0.202(2) 


3.59(7) 


5.8(2) 


Improved 


0.650 




0.384(2) 


3.74(5) 


5.7(1) 


(Landau) 


0.795 




0.304(1) 


3.58(7) 


5.8(13) 




0.905 




0.244(1) 


3.72(8) 


5.7(5) 


Wilson 


1.950 




0.393(2) 


3.17(7) 


5.0(1) 




2.000 




0.355(2) 


3.10(5) 


5.1(1) 




2.140 




0.308(2) 


2.96(5) 


5.4(1) 




2.243 




0.232(2) 


3.05(6) 


5.6(1) 




2.300 




0.203(1) 


3.14(7) 


5.8(5) 




2.400 




0.128(2) 


3.52(11) 


5.8(2) 



TABLE IV. Scalar and tensor glueball masses from anisotropic lattices. 
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FIGURES 



FIG. 1. Static quark potentials on isotropic lattices: (a) plaquette-improved action with 
j3p = 0.730; (b) mean link-improved action with = 0.550; and (c) Wilson action with 
(3w = 1-700. Linear plus Coulomb fits to the on-axis potentials are shown as the solid lines. 

FIG. 2. Effective potential at R = a on the isotropic tadpole-improved lattice with (3p = 0.730. 

FIG. 3. Effective mass plot for the scalar glueball on the isotropic tadpole-improved lattice 
with (3 P = 0.935. 

FIG. 4. Static quark potentials on anisotropic lattices, computed from Wilson loops W x t (•) 
and W xy (□). The potentials in lattice units obtained from W x t have been rescaled by the 
input anisotropy, and the potentials have been shifted by additive constants in order to set 
a s V(R = a s ) = 1. Linear plus Coulomb fits to the Wxt on-axis data are shown as the solid 
lines. Results are shown for (a) the plaquette-improved action at (3p = 0.848; (b) the mean 
link-improved action at (3l = 0.650; (c) the Wilson action, at fiw = 1.950. 

FIG. 5. Effective mass plots for (a) scalar and (b) tensor glueballs on the anisotropic lattice 
with mean link improvement at (3l = 0.795. The result of a single exponential fit to the scalar 
glueball is shown as the solid line, with dotted lines showing the estimated error. 

FIG. 6. Scalar glueball mass versus spatial lattice spacing squared. Data is shown for: the 
mean link- improved action (■); the plaquette-improved action (•); and from the Wilson action 
simulations done here (o). The result of a simulation with an adhoc change to the coefficients 
of the rectangle terms in the improved action, described in the text, is also shown (A). We also 
include results from Wilson action simulations reported in Ref . [^] ( x ) . The star shows the results 
of an a — > extrapolation of published Wilson action data [j3C|l. 

FIG. 7. Tensor glueball mass versus lattice spacing squared. The plotting symbols are the same 
as in Fig. ^. 
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